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DISCOVr.KTKG  HTODIIN  TOTAIXY  LEONTIKF  SUBSTITUTION  SYSTEMS 

by 

George  B.  Dantzlg 

A Totally  l.eontlef  Square  System  Lti  = b is  defined  by  the 
following  properties; 

(a)  Every  column  of  I.  has  at  most  one  positive  element;  b > 0; 

L is  m X m. 

(b)  Lu  = b has  a solution  u ^ 0. 

It  Is  easy  to  see  that  every  row  of  L has  at  least  one  positive 
element  — itu'.eed  there  ate  exactly  m positive  elements,  each  in  a 
different  column.  Other  properties  are:  l is  non-singular,  L ^ > 0 
and  l.u  - b,  ■nL  = c yields  u > 0,  ti  > 0 for  all  h '*  c > 0.  See  (4). 

A.  Totally  Leontlef  Substitution  System  Ax  = defined  by  the 

following  proiv'rt ies : 

(c)  Every  column  Af*,j)  of  A has  at  most  one  positive  component;  b ^ 0 
A Is  m X n . 

(d)  The  linear  program;  Find  Max  z: 

(1)  Ax  = b,  x>0,  ex=*x. 

where  e = (1,  1,  ...,  1),  is  feasible  and  Max  r is  finite. 

Arthur  F.  Veinott,  Jr.  has  made  several  suggestions  that  are  incorporated 
into  the  algorithm  and  proofs  presented,  see  (7). 
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Other  properties  .ire  that  every  column  of  A has  exactly  one 
positive  element;  every  row  of  A has  one  or  more  positive  elements; 
every  subset  of  m columns  with  a positive  elemi^nt  in  each  row  forms 
a feasible  basis.  All  basic  feasible  solut'ons  ari  non-degenerate. 

Main  references  arc  [3],  |7|,  (81,  and  [9j. 

Subst ' ‘-Ct  ion  Classes : Associated  with  any  feasible  basis  B 

to  a general  ’-.neat  program  is  the  set  of  column  ini'ices  j such  that 

A(*,l)  if  si'l'st ituted  for  B(*,i)  forms  a feasible  basis.  is 

non-empty  for  each  1 since  it  includes  j:  A(*,j)  = B(*,i).  In 
general  depends  on  B.  For  the  Leontief  substitution  case,  however, 

it  is  independent  of  the  choice  of  B,  since  the  substitution  class 
clearly  consists  of  all  j that  have  a positive  element  in  the 
same  row.  It  is  convenient  to  have  correspond  to  all  columns  j 

whose  positive  element  is  in  row  i.  We  shall  refer  to  the  latter  as 
"row-standard’ zed"  substitution  classes. 

A system 


(1)  Ax  = b 

where  A is  m n is  a "liidden”  totally  Leontief  substitution  system 
if  there  exists  a nonsingular  matrix  R such  that 

(2)  (RA)x  = (Rb) 

is  a totally  I.eontief  substitution  system. 
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The  Problem.  Given  a linear  system  Ax  = b find  a nonslngular  matrix 
R if  It  exists  such  that  RAx  * Rb  Is  a totally  l.eontief  substitution 
system. 

Constructive  Solution . 

Step  1;  Solve  the  linear  program:  Find  Max  z 

(2)  Ax  = b,  x^O,  ex=z  (Max)  . 

Terminate  if  infeasible  or  unbounded,  or  If  the  optimal  iiasic  fo.aslble 

A, 

solution  is  degenerate.  Let  B be  the  optimal  feasible  basis. 

Step  2:  Determine  substitution  classes  S,,  S. , ...,  S with  respect 

* J z m 

to  B.  Terminate  if  there  is  ambiguity  about  whicii  substitution  class 
a column  belongs  to. 

Step  3:  For  each  i,  for  i = (1,  ...,  m)  determine  multipliers 

R(i,  *)  such  that 

I 

if  j e s^  : 

if  j ? s^ 

i 

Terminate  if  for  any  i = (1,  ...,  m)  the  above  linear  Inequality  system 
is  infeasible.  If  not  terminate,  the  matrix  R whose  rows  are  R(i,  *) 
is  the  solution. 

Proof  Step  1:  If  the  original  system  (2)  is  infeasible  or  unbounded, 

i.e.,  an  optimal  solution  of  the  linear  program  (2)  fails  to  exist, 
then  the  same  is  true  for  every  linear  transform  using  non-singular  R, 
contradicting  the  existence  of  an  R for  (2)  such  that  (d)  holds. 
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R(i,*)  A(*,j)  > 1 , 
< 0 , 

R(i,*)b  > 1 . 


Proof  Step  2:  The  only  way  tht^re  can  be  difficulty  In  determining 

which  columns  an  entering  column  replaces  In  the  basis  Is  the  unbounded 
case  (ruled  out  on  Step  1)  or  there  is  a tie  in  the  ratio  test  for 
selecting  a pivot  in  the  simplex  method.  The  latter  implies  that  the 
feasible  basic  feasible  solution  obtained  after  the  substitution  is 
degenerate.  This  is  ruled  out  since  all  basic  feasible  solutions  are 
non-degenerate  in  a Leontief  substitution  system. 


Proof  Step  3:  Note  that  feasibility  and  the  adjacency  properties  of 

feasible  bases  that  differ  by  one  column  are  invariant  under  non-singular 

linear  transforms.  After  transformation  by  R(i,  *) , assuming  R 

exists,  row  i of  RA  should  have  + components  corresponding  to 

substitution  class  S,  . We  can  permute  the  rows  of  R so  that  S. 

k 1 

corresponds  to  the  "row  standardized"  substitution  class  associated  with 
i of  RA,  namely  R(i,  *)  A(*,  j)  > 0 for  j ^ and  R(i,  *) 

A(*,  j)  0 otherwise;  we  also  have  R(i,  *)*b  > 0 for  i = (1,  ...,  m) . 
Step  III  can  be  solved  by  applying  Phase  I of  the  simplex  method  to  a 
linear  program. 
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